The variation of the electromagnetic coupling and quintessence 



M. C. BentcQ 

Centro de Fisica Teorica de Particulas, Instituto Superior Tecnico 
Avenida Rovisco Pais, 1049-001 Lisboa, Portugal 

R. Gonzalez FelipeQ 
Instituto Superior de Engenharia de Lisboa 
Rua Conselheiro Emidio Navarro, 
1959-007 Lisboa, Portugal and 
Centro de Fisica Teorica de Particulas, Instituto Superior Tecnico 
Avenida Rovisco Pais, 1049-001 Lisboa, Portugal 
(Dated: December 28, 2008) 

The properties of quintessence are examined through the study of the variation of the electromag- 
netic coupling. We consider two simple quintessence models with a modified exponential potential 
and study the parameter space constraints derived from the existing observational bounds on the 
variation of the fine structure constant and the most recent Wilkinson Microwave Anisotropy Probe 
observations. 
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I. INTRODUCTION 

Over the last decade, the temporal and spatial varia- 
tion of fundamental constants has become a very popu- 
lar subject. The motivation partially comes from theo- 
ries unifying gravity and other interactions, which sug- 
gest that fundamental constants could have indeed varied 
during the evolution of the universe It is therefore 
particularly relevant to search for these variations and 
try to establish correlations, if any, with other striking 
properties of the universe, as for instance with dark en- 
ergy. 

From the observational point of view, the time varia- 
tion of the fine structure constant a has been widely dis- 
cussed in several contexts. In particular, from the spectra 
of quasars (QSO), one obtains from the KecK/HIRES in- 
strument [2| 

Aa 

= (-0.57 ± 0.11) x 10" 5 , for 0.2 < z < 4.2 , (1) 

a 

while the Ultraviolet and Visual Echelle Spectrograph 
(UVES) instrument [J 0| implies 



Aa 



(-0.64 ± 0.36) x 10~ 5 , for 0.4 < z < 2.3 . (2) 



The Oklo natural reactor also provides a bound, 

-0.9 x 10" 7 < — < 1.2 x 10" 7 , for z < 0.14 , (3) 

a 

at 95% C.L. @, i, 0. Furthermore, estimates of the 
age of iron meteorites, corresponding to z ~ 0.45, when 
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combined with the measurement of the Os/Re ratio re- 
sulting from the radioactive decay 187 Re — > 187 0s, yield 



Aa 

a 



;±8) x 10" 



at la, and 



at 2cr. 



-24 x 10"' < < 8 x 10" 

a 



(4) 



(5) 



Big bang nucleosynthesis (BBN) also places bounds on 
the variation of a [111 ]: 

-0.05 < — < 0.01 , for 10 9 < z < 10 10 . (6) 

a 

Finally, the 5-year data from the Wilkinson Microwave 
Anisotropy Probe (WMAP) with Hubble Space Tele- 
scope (HST) prior provides the following bound at 95% 
C.L. [l3: 

Aa 

-0.028 < < 0.026 , for z ~ 10 3 . (7) 

a 

Without HST prior this bound is less restrictive 



-0.050 < — < 0.042 , 
a 



(8) 
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at 95% C.L. 

On the other hand, recent observations of high redshift 
type la supernova and, more indirectly, of the CMB and 
galaxy clustering, indicate that the universe is undergo- 
ing a period of accelerated expansion [ 13| . This suggests 
that the universe is dominated by a form of energy den- 
sity with negative pressure (dark energy). An obvious 
candidate for dark energy could be an uncanceled cosmo- 
logical constant [l4| , which however would require an ex- 
tremely high fine-tuning. Quintessence-type models fl5| 
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with one [l6[ or two [TtJ scalar fields, k-essence fl8l and 
the Chaplygin gas with an exotic equation of state [19|,l20( 
are among other possibilities. 

In most of the theoretical approaches, scalar fields are 
present in the theory. Thus, one could expect that the 
two observational facts, namely the variation of a and the 
recent acceleration of the universe, are somehow related. 
Indeed, the coupling of such a scalar field to electromag- 
netism would lead to a variation of the fine structure 
constant HJ. In several contexts [H, 0, the 
above question has already been addressed. 

In this work we shall consider two simple quintessence 
models with a modified exponential potential, proposed 
by Albrecht and Skordis [2g, [2?J ■ Our aim is to restrict 
the parameter space of these models by using the obser- 
vational constraints on the variation of the fine structure 
constant and by imposing consistency with the 5-year 
data from WMAP. 



II. THE MODELS 
A. Quintessence 

The authors of Refs. [26|, [2?J add a polynomial pref- 
actor in front of the exponential potential in order to 
introduce a local minimum in the potential V(<fi), such 
that the scalar field </> gets trapped into it. The potential 
can be written in the form 



■-41 



-X<j> 



(9) 



Since the effect of trapping is equivalent to a cosmologi- 
cal constant, an era of accelerated expansion is eventually 
reached. Notice that this model has already been ana- 
lyzed in the context of the variation of the fine structure 
constant [28|. However, the authors consider a particular 
set of values for the parameters in the potential, whereas 
here we vary these parameters and try to find constraints 
on them. Recently, this potential has also been investi- 
gated in the braneworld context as a model of quintessen- 
tial inflation [2i|. We shall refer to this potential as the 
AS1 model. 

Regardless of initial conditions, trapping occurs in the 
very early universe, when the field enters the attractor 
regime. The tracking properties of the AS1 potential © 
are similar to those of a pure exponential potential. Nev- 
ertheless, due to the presence of the polynomial factor, 
quintessence dominates near the present epoch. Further- 
more, this potential can lead to both permanent and 
transient acceleration regimes. When AX 2 < 1 and the 
field is trapped in the local minimum of the potential, 
permanent acceleration is achieved. On the other hand, 
if AX 2 < 1 and the field arrives at the minimum with a 
kinetic energy sufficiently high to roll over the barrier, a 
transient acceleration regime is obtained 30]. Transient 
vacuum domination also arises for AX 2 > 1, when the 
potential loses its local minimum [30] . The existence of a 



transient regime is interesting from the theoretical view- 
point of string theory, since it can avoid the difficulties 
which typically arise in the S-matrix construction at the 
asymptotic future in a de-Sitter space [3ll. l32l. [33l] . 

A second potential (hereafter referred as AS2 model) 
which can lead to the desired accelerated expansion has 
the form (H, 



V{4>) 



C 



■-fo) 2 + A 



D 



-\<i> 



(10) 



The motivation of such a potential comes from brane 
studies, which indicate that it could arise as a Yukawa- 
like interaction between branes [34J. The behavior of 
the AS2 potential is similar to the AS1 potential during 
the radiation-dominated era. Nevertheless, during mat- 
ter domination it retains much more density than the 
AS1 potential. This is due to the fact that the AS2 po- 
tential has a smoother minimum and a sharper maximum 
when compared to the AS1 one. Furthermore, accelera- 
tion is achieved earlier for the latter potential [27| . 

The evolution of the scalar field is described by the 
equation of motion 



4> + 3Hc/) + V , ((t)) = , 



(11) 



where the dot and prime denote derivatives with respect 
to time and <fi, respectively; 



2 1 

H = ~P, P = Pr + Pm + , 



(12) 



Here p r and p m are the radiation and matter energy den- 
sities, respectively, and p^ — <p 2 /2 + V ((/)). In order to 
integrate equations (jlip and (| 1 2 [1 . it is more convenient 
to rewrite them in the form 

= -3x + V 2 + \x [2x 2 + 7 (1 - x 2 ~ y 2 )] , 

= -a^-x y + \ y [2x 2 + 7 (1 - x 2 - y 2 )}, (13) 

where 7 = wb + 1, tog is the equation-of-state parameter 
for the background (wb = 0,1/3 for matter and radi- 
ation, respectively); J\f = In a, with a being the scale 
factor, and 



f2p 



y 



VP 



"V- ^~T7- (I 4 ) 



XL 
v 



The above equations are integrated from the Planck 
epoch (a ~ 10~ 30 ) to the present epoch (a = a = 1). 
For the fraction of radiation at present we use the cen- 
tral value h 2 — 4.3 x 10~ 5 , which assumes the ad- 
dition of three neutrino species. The matter content is 
such that the correct fraction is reproduced at present, 
Q° m h 2 = 0.1369 ± 0.0037 [35[. Furthermore, we assume 
that early in the radiation era the scalar field is in its 
tracking regime, in which case [36| 



3V, ft ~ — 



(15) 
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B. Coupling with electromagnetism 



Due to the universality of the gravitational interac- 
tions, non-renormalizable couplings of the quintessence 
field <f> to the standard model fields are expected below 
the Planck scale. Following Bekenstein's proposal [2lj . 
one can consider the interaction between the scalar and 
electromagnetic fields in the form 



1 



(16) 



where / is an arbitrary function. Since the variation of 
the electromagnetic coupling is small, one can expand 
this function up to first order in powers of <fi and write 



m = - [i 

a 



■ tw- 



in) 



where £ is a constant. It follows then that the fine struc- 
ture parameter a is given by 



a = [/(£)] ' = ^[1-0 
and its variation is 



Aa 

=0 

a 



'o ) 



(18) 



(19) 



Moreover, the rate of variation of a at present is given 
by 



a 
a 



d<t> 



(20) 



where Ho = 100 h km s _1 Mpc -1 is the Hubble constant. 

At this point, it is worth noticing that any variation of 
a should comply with measurements on the violation of 
the equivalence principle, which translate into the follow- 
ing upper bound on the coupling between quintessence 
and the electromagnetic field [371 ] : 



f < 7 x 10" 



(21) 



III. PARAMETER SPACE CONSTRAINTS 

In our study, we perform a random analysis on the 
parameters of the quintessence potentials described by 
Eqs. © and (flU|) . For the AS1 model we consider the 
two possible late time behaviors: permanent or transient 
acceleration. A viable model should comply with sev- 
eral observational bounds during the different stages in 
the evolution of the universe. A stringent bound comes 
from the amount of dark energy during nucleosynthesis, 
QBBN( Z _ 10 i0) < o 09 com i ng f rom the primordial 

abundances of 4 He and D, which sets the lower bound 
A > 6.7. We scan the full parameter space imposing the 
following conservative bounds at the present epoch: 



0.6 < h < 0.8, 0.6 < Q$ < 0. 
wl < -0.8, q < 0, 



(22) 
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FIG. 1: The minimum value, ^rnin , allowed for the coupling £ 
between quintessence and the electromagnetic field as a func- 
tion of A, for the AS1 model in the permanent acceleration 
case. The horizontal dot-dashed line corresponds to the equiv- 
alence principle bound (pTj) . 
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FIG. 2: Parameter space consistent with all the observa- 
tional bounds, for the AS1 model in the permanent accel- 
eration case. The horizontal and vertical lines correspond to 
the lower bounds on AX 2 and 4>oX, respectively, coming from 
the QSO bounds on the variation of a and the requirement 
|£™„| < 7x ltr 4 . 



where q = —a/(aH 2 ) is the deceleration parameter. 

Our results for the AS1 model are presented in Fig- 
ures Q] to [U The parameter space consistent with all 
the observational bounds is shown in Figs. Q] and [5] for 
the permanent acceleration case. We remark that the 
variation of the fine structure constant a coming from 
the QSO bounds [cf. Eqs. ([1]) and ^j] necessarily re- 
quires a non-vanishing value of the coupling parameter 
£. This imposes a lower bound on this parameter, £ m ;„, 
for the model to be consistent with the observational 
data. Clearly, in order to have a consistent solution, the 
value of |£ m m| should always be lower than the upper 
bound established by the equivalence principle constraint 



4 



- Equivalence principle bound 
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FIG. 3: As for Figure [T] but for the transient acceleration 
case. 
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FIG. 4: As for Figure [2] but for the transient acceleration 
case. 



of Eq. (|2"Tj) . From the figures we conclude that the param- 
eter space of the AS1 model with permanent acceleration 
is restricted to the ranges 



6.6 < A < 21.1 , 
0.01 < AX 2 < 0.97, 
268.1 < O A < 272.3. 



(23) 



We recall that for the AS1 potential the combination 
of parameters AX 2 is useful in distinguishing between the 
permanent and transient regimes, since it determines the 
presence or absence of the minimum in the potential. 
Furthermore, the combination </>oA determines the posi- 
tion of the minimum/maximum or inflection point of the 
potential, tf>± = (1 + 4>qX ± Vl — ^4A 2 )/A, which in turn 
is related to the exit from the tracking regime and to the 
scalar held energy density domination at present. 

The results for the transient acceleration case are pre- 
sented in Figs. [3] and 2J In this case, the equivalence 
principle and QSO bounds on the variation of the elec- 



tromagnetic coupling do not yield an upper bound on 
the parameter A more restrictive than the one already 
imposed by other observational constraints (see Fig. [3]). 
We hnd the following allowed ranges for the parameters 
of the AS1 model with transient acceleration: 



6.5 < A < 21.3, 
0.98 < AX 2 < 1.13, 
270.1 < O A < 272.9. 



(24) 



Our results for the AS2 model are shown in Figures [S] 
and [HI For the sake of comparison with the AS1 model, 
we use the same combination of parameters AX 2 and 4>oX. 
In this case, we find 



6.8 < A < 50.1, 
261.7 < O A < 292.2, 



6.1 x 10" 



< AX < 7.4 x 10" 



(25) 



1.7 x 10~ 14 < C<9.4x 10 2 , 
1.1 x 10~ 14 < D < 8.7 x 10 2 . 



Equivalence principle bound 
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FIG. 5: C, m in as a function of A for the AS2 model. 

Finally, in Figure [7] we give an example of how, for 
particular values of the potential parameters, the AS1 
model fits the observational bounds on the variation of 
the fine structure constant. 



IV. CONCLUSION 

The time variation of fundamental constants is a com- 
mon prediction of theories that attempt to unify the four 
fundamental interactions. The experimental bounds ob- 
tained on the variation of these constants are therefore 
a useful tool for testing the validity of these theories. In 
this work we have studied the implications of the coupling 
of electromagnetism to quintessence fields. In our anal- 
ysis we have considered two simple quintessence models 
with a modified exponential potential (26l . |27] | , with the 
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FIG. 6: Parameter space consistent with all the observational bounds considered, for the AS2 model. 



aim to constrain the parameter space using the observa- 
tional data on the variation of the fine structure constant 
a in combination with the 5-year data from WMAP. Our 
results are summarized in Eqs. (j2"5j) - (j2"5)) . 

We should point out that, in constraining the models, 
the observational data on the variation of a, obtained 
from the quasar absorption systems [cf. Eqs.([T]) and @], 
turned out to be crucial when establishing limits on the 
parameters of the AS1 and AS2 models. This is due 
to the fact that the QSO bounds imply a non-vanishing 
value of Aa/a at redshift 0.2 < z < 4.2, which in turn 
requires a non-vanishing coupling £ between the scalar 



and electromagnetic fields. On the other hand, it was 
not possible to use in our approach the Oklo, meteorite, 
BBN and CMB bounds to put any tight constraint on 
the model parameters, since these observational bounds 
are all consistent with no variation of a at their corre- 
sponding redshift scales. 
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